fig = figure;
plot(xi_oecd,f_oecd,'- k','LineWidth',2);
hold on;
   plot(xi_nonoecd,f_nonoecd,'-- k','LineWidth',2);
hold off;
lgd=legend('OECD','Non-OECD','Location','north');
% lgd.FontSize = 30;
% ax = gca;
% ax.FontSize = ax_fs;
xlim([0.5 1.8]);
xlabel('{\sigma}_c');
orient landscape;
print(fig, [figdir, 'Figure_3_d' mod_name], '-dpdf','-bestfit');
close(gcf);